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INTRODUCTION 

The  recovery  of  temperature  profiles  of  the  earth's  atmosphere  by 
a  top-side  sounder  in  the  15  pm  region  has  been  investigated.  A  modeled 
radiative  transfer  of  COg  through  a  structured  atmosphere  is  created  and 
analyzed  by  a  least  squares  method  in  the  presence  of  varying  amounts  of 
random  noise.  This  is  done  both  in  the  spectral  domain  using  either 
a  135  cm"^  window  centered  at  15  pm  or  by  using  10  selected  spectral 
channels  0.5  cm“^  wide  or  in  the  interferogram  domain  using  either  the 
entire  continuous  region  from  zero  path  to  an  optical  path  difference 
corresponding  to  the  spectral  domain  resolution,  or  by  using  selected 
portions  (partial  interferograms)  to  recover  temperature  profiles.  The 
recovered  profiles  and  their  root  mean  square  errors  are  then  compared  both 
in  the  spectral  and  interferogram  domains  in  order  to  evaluate  the 
apparent  advantage  in  using  partial  interferograms  as  first  proDosed  by 
Kyle  (1977)  and  investigated  by  Smith  et  al.,  (1979)  and  at  Santa  Barbara 
Research  Center  (1981)  and  Perkin- Elmer  (1982). 

First  a  computation  of  synthetic  line  spectra  for  the  nine  strongest 
bands  of  CO^  at  15  pm  was  done.  These  are  general  spectra  parameterized 
by  temperature,  pressure,  column  density,  and  general  line  widths.  The 
model  then  convolves  a  Lorentz,  or  a  lorentz-Gaussian  (Voigt)  profile 
with  the  line  spectra  and  forms  the  transmission  and  absorption  as  a 
function  of  wavenumber. 

The  radiative  transfer  for  multiple  atmospheric  layers  (up  to  14 

distinct  layers)  was  then  computed  and,  for  spectral  domain  studies, 

convolved  with  an  instrument  function,  either  triangular,  sine,  or  preferably 
? 

sine  .  Fourier  transforms  of  the  separate  and  composite  atmospheric  layer 


spectra  were  also  done,  since  they  are  required  for  the  interferometer 


domain  analysis.  The  final  result  is  an  artificial  spectrum  representative 
of  that  observed  by  a  satellite  or  balloon-borne  interferometer  looking  down 
through  a  radiating  atmosphere  to  a  hot  earth  surface.  The  possible  gain 


advantages  are  then  discussed  and  conclusions  drawn. 


Radiative  Transfer  in  a  Plane  Parallel  Atmosphere: 


We  assume  a  plane  parallel  atmosphere  with  no  scattering,  and  in 


local  thermodynamic  equilibrium,  so  that  a  particular  layer  may  be 


characterized  by  a  pressure,  temperature,  and  density,  and  a  constant 


mixing  ratio  of  COg*  Under  such  a  condition,  the  radiative  transfer 
calculations  are  done  using  the  expression  for  the  intensity  of 
radiation  I(o)  emerging  from  the  top  of  the  last  layer  (looking  at  nadir) 


Aop 

=  l  B(o , 


bottom 


T)d-r(a)  +  B(a,Tg)xg(o) 


where  x(o)  is  the  wavelength  dependent  transmission  from  a  given  layer 
at  height  z  to  the  top  layer,  and  B(o,t)  is  the  black  body  function  for 


a  given  layer  (z  dependent),  and  the  subscript  g  is  the  ground  contribution. 


Of  course  our  integration  is  done  with  only  a  finite  number  of  layers 


(up  to  14  homogeneous  layers  of  about  2  km  thickness  plus  the  ground 
contribution)  so  the  expression  becomes  a  summation,  so  that  we  have 


I(o)  “  22  “iB(a»TiMTi,i  _Ti,i+i)  +  B(°Jghu 

1=1 


where  are  numerical  weights  for  integration. 

Because  of  the  nature  of  the  transmission  function  in  the  form  of 
an  exponential  whose  argument  is  proportional  to  the  thickness  and 
optical  absorption  coefficient  for  a  particular  layer,  the  total 
transmission  from  level  1  to  level  i  is  the  product  of  the  individual 
transmissions  within  each  homogeneous  layer,  so  we  have,  for  example, 

t13  —t j2t23*  t14  =  t12t23t34  =  T13T34’  etc*  In  this  way  we  build  up 
the  total  expression  for  1(a)  for  a  given  temperature  profile.  Such  a 

calculation  is  illustrated  from  Figures  1-11  for  ten  atmospheric  levels 

and  a  ground  contribution,  and  represents  the  contribution  of  each  layer 

to  the  radiance  at  height  30  km  for  a  down-looking  instrument.  Figure 

12  shows  the  total  radiance. 


Least  Squares  Procedure  for  Recovering  Temperature  Profiles. 


ectral  Space. 

In  order  to  derive  a  temperature  profile  from  experimental  measure¬ 
ments  of  intensity,  the  most  direct  method  is  to  assume  a  reasonable 
"most  probable"  profile  appropriate  for  the  season  and  region  and  work 
with  the  difference  spectra,  with  deviations  indicating  actual  temperature 
variations  from  the  assumed  standard  profile.  For  small  changes  in 
temperature  a  Taylor  series  expansion  is  sufficient,  and  we  retain  only 
the  first  differences.  Thus  an  intensity  profile  giving  a  nadir  intensity 
of  1(a)  may  be  represented  as 


I(o)  .  I0M  ♦  £  [%Il|  4Tj 

<■'  L  >1 
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where 


ai(a.T) 

3T 


■f(aJ-Ti)  S 


fji 


la=a. 

J 

T=Ti 


(4) 


We  minimize  ^(AT^)  in  the  least  squares  sense  as  (for  a  set  of  discrete 
wavenumbers,  a.,  l<j<M) 

J  ” 


(5) 


The  set  AT^  is  chosen  such  that  4>  is  a  minimum,  d<f>=0,  or  since  the  set 
ATi  is  independent 


H  r 

3*  -  0  =  (-2)  £ 

W  j=i 


N 

(l(°j)”I0(aj))  '  2  f(aj’Ti)ATi 

i=l 


f(°j  *Tk)  -(6) 


Letting^Oj)  ■  I(oj)  -  I0(Oj)  *  the  difference  between  the  observed 
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spectra  and  the  spectra  from  the  standard  profile,  we  have 


£iiaj)f(VV  ■ 

j=l 


H 

jal 


AT, 


(7) 


This  set  comprises  N  equations  for  the  profile  differences  from  a  set  of 
M  points  of  spectra  data. 

We  define  the  following  matrices,  first  the  MxN  matrix  of  coefficients. 


fll  f12  ***  flN 

f21  f22  *•' 

•  • 

•  • 

•  • 
fMl  *  *  *  fMN 


(8) 


Next  the  temperature  difference  column  Nxl  matrix. 


’iTi" 

AT  ■ 

AT2 

and  the  intensity  difference  matrix, 

• 

yy 

-y 

1 

Then  the  expression  for  the  temperature  differences  can  be  written  in 
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matrix  form  as 


fT/  =  fTf  AT  ,  (10) 

where  is  the  transpose  operation  on  f,  and  makes  the  product  matrix 
f*f  square,  and  possessing  an  inverse. 

By  including  a  "noise"  term  in  the  expression  (  3 )  we  can  determine 
the  variance  in  the  temperature.  Thus,  adding  a  noise  term  n,  we  have, 
in  matrix  notation 

J  +  n  =  f AT  ,  (11) 

where  n  is  an  Mxl  column  matrix.  Then  <p,  the  quandratic  expression  to 
be  minimized,  may  be  written  as 

y 

<p  =  [/+  n  -  fAT]  s[j+  n  -  fAT]  ,  (12) 

where  S  is  a  symmetric  weighting  matrix  to  be  determined. 

Multiply  (11)  on  the  left  by  fTS, 

fTsf  +  Vsn  =  flSf  AT  •  (13) 

Now  fTSf  usually  has  an  inverse,  called  the  Moore-Penrose  pseudo- Inverse 
with  weight  S.  Thus  we  may  solve  for  AT  as 

AT  «  (fTSf)_1fTsy  ♦  (fTSf)_1fTSn  ’  (14) 

The  "best  estimate"  AT  of  AT  is  the  first  term. 


(15) 


AT  = (fTSf)-1  fTS j 

and  the  "estimation  of  the  error"  is  then 


(AT  -  AT)  =  (fTSf)-1fTSn 


(16) 


Now  the  nature  of  n  is  stochastic,  and  the  only  meaningful  measurements 
of  our  error  are  its  expectation  in  the  probability  sense.  So  if  we 
define  a  functional  operation  as  the  expectation  of  a  function  g,  as  E(g), 
and  use  a  probability  density  function  p,  such  that 


and  where 


p(v)dv  =  1 
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p(v)dv 


» 


is  the  probability  of  an  event  occurring  between  a  and  b,  then 


(17) 


(18) 


oo 

g(v)p(v)dv  (19) 

CD 

is  the  "expectation  value  of  g".  Then  the  variance  of  the  estimation 
error  may  be  expressed  as 


14 


E(AT  -  AT) (AT  -  AT)T 


=  Ei 


|(fTSf)‘1fTSn] |(fTSf)-1fTSn]  1 


*  E^  (fTSf)"1fTSn  nTSf(fTSf) 


(fVrVs  E(nnT)  Sf^Sf)’1 


(20) 


Now  E(nn  )  is  the  variance  of  the  matrix  Q,  the  noise,  and  we  may  minimize 
<j>  and  e[(AT  -  AT) (AT  -  Sf)^]  simultaneously  by  choosing  our  weighting 
matrix  to  be  the  inverse  of  the  noise  matrix,  Q.  Thus, 


S  =  Q"1 


(21) 


Then 


E  [  (AT  -  AT) (AT  -  A?)1]  =  (fW)'1 


(22) 


and  the  "best  estimate"  of  the  temperature  profile  is 


AT 


=  (fTQ-1f)-1 


(23) 


Now  in  general  the  noise  matrix  nn  is 

n-, 


'I 


n, 

t 

n. 


(njn2n3...nu) 


n_T  _ 
nn  * 


r  nj  njn2.. 
t^i  n2  n2^ 


(24) 
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For  uncorrelated  noise,  the  expectation  value  of  nn  leaves  only  the 
diagonal  elements,  as 


Q  =  E(n  n  )  * 


E{nj)  0  0  . 

0  E(n?2) 


We  may  further  reasonably  assume  the  noise  to  be  distributed  uni  formally, 
so  the  diagonal  terms  all  have  the  same  value,  so  Q  becomes  proportional 


to  the  unit  matrix. 


1  0  0 


Q  =  °n  I  0  1  0 


where  a  is  the  variance  of  the  noise.  This  being  the  case  we  can  then 
n 

write  our  expressions  for  the  estimation  of  the  temperature  variance  as 


E  |^(AT  -  AT) (AT  -  AT)T]  =  o2(fTf)_1  , 


and  the  best  estimate  of  the  profile  difference  as 


AT  =  (fTf)_1  fT  f 
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The  Recovery  of  Temperature  Profiles  from  Interferograros. 

The  Interferogram  E(x)  is  essentially  the  Fourier  cosine  Integral 
transform  of  the  spectrum,  I(o),  as  (Vanasse  et  al.,  1971) 

r* 

E(x)  *  2e  I  I(o)cos  2i rox  do  ,  (29' 

•'o 


where  e  is  the  efficiency  of  the  interferometer. 

In  practice  the  spectrum  is  computed  at  2N  discrete  points,  as  shown 


The  integral  is  then  related  to  the  complex  Fourier  transform  as. 


first 


cos  2troxdo 


(30) 


and  then  converted  to  a  discrete  sum  by  multiplying  I (a)  by  the  Dirac 
comb  function  and  by  band  limiting  by  a  rectangular  function,  as 


f  1(0)  comb 
-no 


Rect 


p-  (gp-r ) 

(2N-1)Ao 


e-i2irox  do 


(31) 


where 


Ao 


<°N  '  °0> 
N 


(32) 


is  the  wave  number  sampling  increment.  Also  note  that  the  nth  sample 
point,  on,  is  expressed  as 


o„  s  on  +  nAo  , 
n  o 

and  that  the  width  of  the  rectangular  function  window  is 


2(<jn_i  -  o_M)  =  (2N-1)ao 
We  may  complete  this  as 


(33) 


(34) 
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E(x)  *  Ao  I  I(o)  Rect 


o-fap-y 

(2N-l)io 


Aov  +oo 


V  6(O-o0-nlo)e-,2"ox  do 


n=-oo 


^i2"x<vn40) 

n=-oo  ^  -* 


Since  the  definition  for  a  rectangular  function  Rect  x  is 


Rect  x 


!*  ’  -7 


0*  otherwise 


then  the  sum  is  limited  to  -N  ^  n  _<  N  -  1  . 

In  addition  the  optical  path  difference,  x,  is  also  measured  at  discrete 
intervals,  so  xm  =  mAx,  m  integer,  so  we  may  write 


N-l 

.  V'  -i2iTmAx(oft+nAo) 

E(mAx)  *  Ao  Ifag  +  nAo)e  u 


-127rAxonm 


5  «..* 


-i27rAxAomn 

nAo)e  •  (37) 
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-i 


b 


For  I(o-  +  nAo)  *  I(oA  -  nAo),  and  AxAo  =—  (making  use  of  the  sampling 
0  0  2N  V 


theorem  Ax  B 


1“) 


we  have 


2(«H-r  °o> 


E(mAx)  *  Acre 


i2irAx0Qm 


-1 


S  Ho0+ 


nAa)e 


-i2TrAxAomn 


n«  -N 


-i2irAxAomn 

+i(c0)  +  2-*  I{ao +  1,40)6 

n  =  1 


=Aa  e 


•i  2irAxagm 


N-l  /  + 

fl(a0-nAo)e 
n=l  \ 


iirmn 

“IT 


+  I(og+nAo)e 


mmn  x 

■) 


+  I(ctq)+  I  (oq+NAo) (-1) 


m 


=  Ao  e 


-i2"irAxagm 


N-l 


I(Og)  +1  (oq+NAo )(-l)m+2  ^I(ag+nAo)cos' 
•  n=l 


Timn 


-i2irAxanm 

In  addition,  by  choosing  Axog  =  h,  integer,  e  =1,  (which  makes 
oQ  an  even  multiple  of  the  bandwidth  (oN_j-  oQ)  or  oQ  *  2h  (c^  -  0Q) 
and  so  we  finally  have  for  the  discrete  Fourier  cosine  transform 


(38a) 
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E(mAx)s  JT  1(c)  = 


N-l 

I(aQ)  +  2  I(  Cq+  nAo)  cos  -^-^+I(oq+Nao)  ( — l)m 


|Aa . 


Profile  Recovery  from  Interferometer  Data: 

By  applying  the  discrete  Fourier  cosine  transform  to  the  spectral 
space  matrix  equation  (11  ),  we  have 

^[/(o)+  n(o)]  [f (a)  ]  AT  .  (39) 

We  retain  the  linear  nature  of  the  analysis,  since  the  two  column  matrices 
and  n(o)  and  the  rectangular  matrix  f(o)  transform  directly,  and 

gives 


£(x)  +  N(x)  =  F(x)  AT  •  (40) 

Then  by  applying  a  similar  analysis  as  before  we  have  for  the  estimation 
of  the  temperature  variance 

E  |  (AT  -  A?) (AT  -  AT)T]=  c^(FTF)-1  ,  (41) 

and  the  best  estimate  of  the  profile  difference 

AT  -  (FTF)-J  FT£ 

2 

where  o^  is  the  variance  of  the  delay-space  noise. 
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(38b) 


(42) 


How  Noise  Transforms; 

We  use  Parseval's  theorem  as 


n2(a)da 


N2(x)dx 


(43) 


where  N(x)  is  the  Fourier  cosine  transform  of  n(cr).  For  the  discrete, 
band  limiting  procedure  previously  applied  to  the  Fourier  cosine  transform, 
we  get 


n2(a0)  +  2  2^  n2(<*o  +  nA°)  +  1)2  (<VNA°) 

n=l 


A  X 
N 


N-l 

N2(0)  +  2  N2(nAx)  +(ag+NAo) 

n=l 


(44) 


These  expressions  are  just  the  noise  variances,  and  so  we  have 


Ao 


A  x  a 


N 


(45) 


This  Is  used  to  compare  the  effect  of  equivalent  noise  in  the  two  transform 
domains  when  we  make  our  numerical  analysis. 
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Some  Numerical  Relations  and  Parameters  Used  in  Our  analysis: 

We  present  here  some  parameters  used,  along  with  pertinent  algebraic 
relations.  The  number  of  points  to  model  the  spectral  domain  Is  N  =  4096, 
The  beginning  wavenumber  Og  =  600  cm  1  and  the  end  wavenumber  o^  j=  750 
giving  a 

Ao  =  ---‘i  —  =  .036621  cm"1  .  (46) 

"  4096 

Since  by  the  sampling  theorem  the  maximum  O.P.D.  delay  L  Is 


a  -  1 

Ao  — —  * 


then  L  =  13.653  cm 


(47) 


The  delay  space  Increment  is 


Ax. 


2{aN-r  °o) 


=  .00333+  cm 


(48) 


The  relation  between  the  noise  variances  may  be  written 


o2  . 
N  " 


N 


2 


(49) 


The  computer  program  returns  the  Fourier  cosine  transform  without  the 
Ao  term,  thus  we  have 


o  o 

So,  for  example,  r.m.s.  noise  of  =  10  watts/cm  -sr  is  a  spectral 


r.m.s.  noise  of 


lO-8 


watts 
~2 


cm  -sr-cm 


-1 


=  3.0170  x  10’®  watts/cm^-sr-cm" * 


3.3146 

and  an  r.m.s.  noise  returned  by  the  computer  of  ( 90.5097 ) ( 3 .01 70x10’®)  = 
2.7307  x  10’7  watts/cm^-sr. 
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Actual  Profile  Recovery  for  Simulated  Spectral  Models. 

Table  3  shows  the  location  for  ten  "optimal"  filters  suggested  by  Kaplan 
et  al.  (1977).  These  were  applied  to  the  spectral  region,  using  each  filter 
region  as  a  sum  over  the  window  to  give  the  equivalent  response  for  ten 
separate  radiometers,  and  then  the  spectral  space  covariance  matrix  (fTf)-1 , 
Eq.  (  27  )  was  found.  The  ten  atmospheric  layers  captioned  in  Figures  1-10 
were  used  for  the  spectral  space  filters.  Table  15  shows  the  covariance  error 
matrix,  with  the  diagonal  elements  being  the  coefficients  for  the  temperature 
variance  to  be  multiplied  by  the  variance  of  the  noise  present  from  detectors, 
etc.  The  r.m.s.  noise  temperature  for  the  ten  levels  becomes 

AT  =  (1.95  x  1 08 ) ( 3 .017  x  10"9)  =  .587  °K  (53) 

_g  wflt 

where  we  assume  that  3.017  x  10  --—a"-  r  r.m.s.  noise  is  present. 

cm  -sr-cm~L 

A  second  analysis  uses  the  full  spectral  information  from  600-750  cm”1, 
assuming  fourteen  atmospheric  layers.  The  noise  temperature  error  matrix 
(fTf)_1,  is  shown  in  Table  16  with  an  r.m.s.  noise  temperature  of  .15°K,  again 
assuming  the  same  amount  of  noise  as  previously. 

For  the  transform  domain  simulation,  Fourier  transforms  of  the  separate 
and  composite  atmospheric  layer  spectra  are  performed  and  the  error  covariance 
matrix  (FTF)-1  is  constructed,  using  either  the  entire  transform  optical  path 
difference  interferogram  (0-13.653  cm),  or  various  partial  interferograms,  as 
displayed  in  Tables  4-14.  Note  that  in  comparing  the  two  domains  we  use  the 
equivalent  amount  of  noise  for  each,  and  are  not  considering  the  obvious 

V'g  -  <7  .  ’ 

,  by  observing  all 

elements  at  the  same  time  (as  with  an  interferometer  versus  a  grating 


spectrometer).  This  is  done  to  emphasize  the  information  content  of  the  two 
domai ns . 


Gain 

Figure  41  shows  the  r.m.s.  noise  temperature  given  by  the  matrix  (FTF)-1 
for  various  partial  interferograms  with  10"8  watts/cm2-sr  noise.  The  most 

promising  region  appears  to  be  the  continuous  i nterferogram  (0-2.167  cm), 
which  has  a  r.m.s.  value  of  .276  °K  for  650  points  of  data.  If  we  were  to 
spend  an  equal  amount  of  time  measuring  the  partial  interferograms  from  either 
(.500  -  .833  cm)  with  an  r.m.s.  noise  temperature  of  1.011  °K  for  100  points, 
or  the  region  (1.167  -  1.500  cm)  with  an  r.m.s.  noise  temperature  of  .969°K 
for  100  points,  we  would  expect  the  r.m.s.  noise  to  decrease  by  the  ratio 
VlOO/650  =  1/2.550  ,  and  thus  we  could  expect  to  determine  the  temperature  for 
the  (1.167  -  1.500  cm)  region  to  within  =  .380°K,  which  is  still  not 

as  good  as  the  .276°K  r.m.s.  for  the  0-2.167  cm  region.  Figure  42  shows  th? 
comparison  for  the  (0-2.167  cm)  region  with  the  (.500  -  .833  cm)  region  using 
equal  measuring  times.  Thus,  under  the  limited  cases  studied  and  subject  to 
the  model  accuracy  we  find  no  gain  advantage  for  using  partial  i nterferograms 
in  the  15  ym  region. 
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Figure  21.  The  function  B'-.goft 
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Figure  24.  The  function 


Figure  27.  The  Fourier  transform  of  Figure  13,  FT  B'290°^T30-2  "  T30>0 
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PER  0 1  V I S I  ON  =  1  X  10- 


Four  ler  transform 
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Figure  31.  The  Fourier  transform  of  Figure  17,  FT  B'242o(T30-10  '  T30-8 
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Fourier  transform  of  Figure  18,  FT  B‘  „(t,„  -833  cm 
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Figure  45. 


Computation  procedure  used  to  generate  transmittance. 
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BAND  CENTER  UPPER  LEVEL  IQWfR  LEYBL.  ISOTOPE 


BAND  STRENGTH 


667.379 

01101 

00001 

626 

82580 

647.058 

11102 

10002 

626 

222 

648.484 

01101 

00001 

636 

860 

688.678 

11101 

10001 

626 

148.8 

622.368 

01101 

00001 

628 

330 

668.109 

03301 

02201 

626 

382.4 

667.750 

02201 

01101 

626 

6488 

618.033 

10002 

01101 

626 

1436 

720.808 

10001 

01101 

626 

1853 

TABLE  2.  STANDARD  ATMOSPHERE  —  MIDLATITUDE  SUMMER 


LAYER  HEIGHT 

TEMPERATURE 

PRESSURE 

ATM  -  CM 

0 

291 

1.000 

0-2 

290 

.792 

58.9 

2-4 

279 

.620 

47.5 

4-6 

267 

.481 

37.6 

6-8 

255 

.367 

30.0 

8-10 

242 

.277 

23.1 

10-12 

229 

.206 

17.8 

12-14 

216 

.151 

13.4 

14-16 

216 

.110 

9.81 

16-18 

216 

.0802 

7.01 

18-20 

217 

.0587 

5.00 

20-22 

219 

.0431 

3.60 

22-24 

222 

.0318 

2.60 

24-27 

224 

.0216 

2.36 

27-30 

234 

.0130 

1.96 

30- 

243 

.0130 

TOTAL  260.6  Atm. 
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TABLE  4.  COVARIANCE  ERROR  MATRIX  (FT)" 
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RMS  (DIAGONAL)  =  0.702371E+06  Window  0.4096,  (0-13.653  cm) 


TABLE  5.  COVARIANCE  ERROR  MATRIX  (FT)' 
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RMS  (DIAGONAL)  =  0.153622E+07  Window  50-650,  (.167-2.167  cm) 
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0,84337E+11  -0.10241E  +  12  0.50030E  +  11  0.96588E 

RMS  (DIAGONAL)  =  0.158840E+07  Window  100-650,  (.333-2.167  cm) 
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